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We study the ordered phases and the phase transitions in the stacked triangular anti- 
ferromagnetic Ising (STAFI) model with strong interplane coupling modeling CsCoCls 
and CsCoBra. We find that there exists an intermediate phase which consists of a single 
phase of so-called partial disordered (PD) type, and confirm the stability of this phase. 
The low temperature phase of this model is so-called two-sublattice ferri magnetic phase. 
The phase transition between the PD phase and two-sublattice ferri magnetic phase is 
of the first order. This sequence of the phases is homomorphic as that in the three- 
dimensional generalized six-state clock model which have the same symmetry of the 
STAFI model. By studying distributions of domain walls in one dimensional chains 
connecting layered triangular lattices, we clarify the nature of the phase transition and 
give an interpretation of little anomaly of the specific heat. 
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1. Introduction 

The hexagonal ABXs-type crystals CsCoCls and CsCoBrs have attracted a great interest be- 
cause of the successive phase transitions due to strong frustration. For a long time, a large number 
of experimental and theoretical studies have been done for this type of substances. The magnetic 
ion Co^^ has fictitious S = 1/2 Ising spin. These substances consist of a quasi-one-dimensional spin 
chains because the exchange interaction along the c-axis is much larger than that in the c-plane,^-* 
and they are coupled via antiferromagnetic interaction in the c-plane, This type of substances is 
expressed by the stacked triangular antiferromagnetic Ising (STAFI) model. 

The neutron scattering experiments on CsCoCls^^ and CsCoBrs^^ suggested that this type of 
substances has the two magnetic phase transitions. In order to explain those experimental results, 
the theoretical study on this model was done by Mekata.^^ Due to the strong interaction along 
the c-axis, one may assume that, in the ordering process, the spins in a chain can be regarded as a 
magnetic moment, and these magnetic moments form a triangular lattice. Mekata studied the phase 
diagram of this two-dimensional model by using a mean field approximation (MFA). He introduced 
the next-nearest neighbor interaction in the c-plane to explain the successive phase transitions 
observed in experiments. He found that this model has three type of ordered phases. That is, 
as the temperature is reduced, there appear the partial disordered (PD) phase, three-sublattice 
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ferrimagnetic (3FR) phase and two-sublattice ferrimagnetic (2FR) phase which are schematically 
depicted in Fig.l. In the PD phase, the spins on one of the three sublattices are disordered and 
the spins on the other two sublattices are ordered antiferromagnetically. In the 3FR phase, the 
spins on two sublattices are ordered in the same direction and the spins on the last sublattice is 
ordered in the opposite direction. The magnitudes of the sublattice magnetizations are mutually 
different. In the 2FR phase, the magnitudes of the parallel sublattice magnetizations are equal. 
The temperature dependence of the magnetic structure factors are in good agreement with the 
experimental observations. However, in the experiments, the specific heat has no anomaly at the 
low temperature transition point. As to the nature of the intermediate phase, many possibilities 
have been proposed, and the nature of the intermediate phase has been mystery for a long time. 

In the numerical studies, it had been difficult to simulate the STAFI model by using a single 
spin flip Monte Carlo (MC) method owing to its frustration and strong one-dimensional anisotropy. 
Recently, Koseki and Matsubara conquered this difficulty by introducing the cluster heat bath 
(CHB) algorithm.^) They obtained the results that agree with the experiments. Their results 
reproduce experimental data about the magnetic structure factors and the specific heat (see Fig. 2). 
However, in their results, there were large fluctuations in the intermediate temperature phase 
corresponding to the PD phase. Thus, the nature of the intermediate phase and property of the 
low temperature phase transition have been not clear. In particular, at low temperature, the 
correlation length along the c-axis in the disordered chains is very long and it seems to diverge in 
simulations in small system sizes. Therefore, it could be possible that the low temperature phase 
transition may not be a transition in equilibrium but a crossover for the finite size effect. 

From the view point of the symmetry, the STAFI model has the Ze symmetry^) and the 
behavior of the phase transition of this model may be explained by the three dimensional six-state 
clock model which has also the Zg symmetry. The simple six-state clock model has been found to 
have only one type of symmetry broken low temperature phase.^^ On the other hand, it turned out 
that a generalized six-state clock (6GCL) model has an intermediate phase where two neighboring 
states are mixed (incompletely ordered phase, lOP). In this paper, we propose that the structure of 
phases in the 6GCL is the same as that of the present STAFI model. In the 6GCL, phase transition 
between the intermediate phase the paramagnetic phase belongs to the universality class of three 
dimensional XY model, and the phase transition between the lOP and the low temperature phase 
is of the first order.^) Therefore, we expect, from the correspondence between the STAFI model 
and the 6GCL model, that the intermediate phase on the STAFI model is the PD phase, and the 
phase transition between the PD phase and the low temperature phase is of the first order. It has 
been known that the intermediate temperature phase of the 6GCL model consists of a single phase 
of lOP, and thus we expect a single intermediate temperature phase, i.e., the PD phase, exists 
in the STAFI model. That is, the complicated phases so far proposed or the 3FR phase do not 
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appear. In the present paper, we investigate characteristics of the low temperature phase transition 
by using the MC simulation and cluster mean field approximation. 



2. Model and Numerical Results 

The Hamiltonian of the STAFI model is given by 



nnn 



i (ix,v) (ix,\) 

where Si^i, = ±1. ^"-^ and run over nearest neighbor pairs and next nearest neighbor pairs 

in the c-plane of the stacked triangular lattice (Fig. 3), respectively. The sign of the exchange 
interactions Jq and Ji are negative while J2 is positive. We set exchange interaction to \Ji/Jq\ = 
3 X 10~^, \ J2/ Jq\ = 1 X 10^^ which are the same as those used by Koseki and Matsubara.^-* We use 
lattices of the size 6n x 6n x 250n with n = A ^ 7 . We impose the periodic boundary condition for 
all directions. We perform MC simulation with the CHB algorithm. The initial 1000 x n steps 
are used for equilibrating process. After this, 5000 x n steps are used to take thermal averages of 
physical quantities. 

We calculate the following quantity:^' 

-1 / 



M, = ^co.(^6Un-'(^^Jj^, (2) 
where 

nix = J^ii - ^(^2 + ma), (3) 

my = ^{mz-m2) (4) 

and rrii (i = 1,2,3) is the sublattice magnetization of i-th. sublattice. The present model can be 
transformed into a six-state clock model as depicted in Fig.4. The quality 

^ = tan~'''(m3;/ma;) (5) 

is the angle in the transformed model. Then, Mg = (cos 6^) detects the order of the six-state clock 
model and it takes 1 for the 2FR phase and —1 for the PD phase. When the transition is of the 
first order, this quantity is scaled by the system size as 

(cos 6^) ^ /(A;b(T - rc)nVJo), (6) 

around the transition point. Here d is the dimension of the system, Tc is the transition temperature, 
and f{x) is a finite size scaling function. We obtained a good scaling plot with the /cbTc/Jq = 0.188 
as depicted Fig.5. Therefore, we estimated that this transition is of the first order and the transition 
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temperature is A;b7c/Jo = 0.188(2). In the high temperature region beyond Tq, Mg goes to — 1, 
which means that the system is stable in the PD phase. The macroscopic fluctuation which was 
proposed in the previous study^^ does not appear. The large fluctuation is attributed to large 
fluctuation due to the apparent XY behavior which is also observed in the 6GCL model.^^ The 
3FR phase is not observed. In the present study, we used a larger sized lattices. Therefore we 
succeeded to detect bulk structures in the large fluctuation. In Fig.5, we see a good convergence 
of date for n = 6 and n = 7. 

Although the low temperature phase transition is of the first order, the specific heat shows no 
anomaly at the transition point as depicted in Fig.6. It means that the latent heat of this transition 
does not exist or at least is very small. The difference between the PD phase and the 2FR phase is 
whether the disordered chains exist, or not. The energy difference between the ordered chain and 
the disordered chain is considered to be given mainly by the number of domain walls rid in a chain. 
Because of the long correlation length along the c-axis, is very small. In order to study how the 
disordered chain changes to the ordered chain, we examine the domain wall density 

Pi = nd,i/L (7) 

in the i-th sublattice. Here L is the length of the chain. We show the temperature dependence 
of the density of domain walls for the case of n = 7 in Fig.6. In the PD phase, the inter-chain 
interactions cancel out each other on the disordered chain. Therefore this chain can be regarded 
as a one dimensional Ising chain with interaction Jq, when we ignore the next nearest neighbor 
interaction J2 because IJ2/J0I = 1 x 10~^. The domain wall density of the disordered chains pi is 
roughly represented by that of the one-dimensional Ising model 

In Fig.6, we find a discrete change around k^T/J^ ~ 0.188. The region k^T/J^ > 0.188 corresponds 
to the PD phase. There, pi is a little smaller than pid which is attributed to the next nearest 
neighbor interaction. The other two chains are ordered chain and, the densities of the domain 
wall, p2 and ps are small, and ptot = Pi + P2 + Ps — Pi- On the other hand, in the 2FR phase 
{k^T/Jo < 0.188), the domain walls move to the two less-ordered chains (pi — P2)- Thus, the 
configurations of the domain wall distribution changes discontinuously, which is consistent with 
the first order phase transition. Here it should be noted that the ptot changes very little at the 
transition point. Thus, the latent heat (or the energy differences) between the PD phase and the 
2FR phase due to the domain wall is very small. Therefore, the dominant contribution to the 
latent heat comes from energy due to the next nearest neighbor interactions. For CsCoCla and 
CsCoBr3, the contribution due to the next nearest neighbor interaction J2 is much smaller than 
that of the domain wall (in the present study \J2/Jo\ = 1 x 10~^). Thus, it is difficult to measure 
of the anomaly of the specific heat in the experimental and numerical observations. 
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3. Cluster Mean Field Approximation 

As mentioned in the introduction, as to the existence of the 3FR phase, there is a discrepancy 
between the result of the MFA and of the Monte Carlo method. In order to resolve this disagree- 
ment, we study the effect of fluctuation in the MFA. In order to take the effect of fluctuation 
into account, we adopt so-called cluster mean field approximation where we investigate clusters of 
spin with the mean field representing the interaction with the outer spins. Therefore, fluctuations 
of spins inside the cluster are taking into account. We takes the cluster depicted in Fig. 7. The 
temperature dependence of the sublattice magnetizations are plotted in Fig. 8. In the model (a) we 
reproduced the Mekata's result. In the model (b) the region of the 3FR phase decreases, and in 
the model (c), it disappears. Therefore, we conclude that the 3FR phase does not stable against 
fluctuation. The suppression of the 3FR phase has been reported in the studies using other ex- 
tended mean field theory.iO'") If the 3FR disappears, the PD phase can not change to the 2FR 
phase smoothly. Therefore we have the first order phase transition. These results are consistent 
with the results of MC simulation in the previous section. 

4. Conclusion and Discussion 

We study properties of ordered phases of the STAFI model. At low temperature, the present 
model has 2FR phase. Between the 2FR and the paramagnetic phase at high temperatures, there is 
an intermediate phase of the spin structure of the PD phase. There is no the 3FR phase predicted 
by a simple MFA. This structure of phases is homomorphic to that of the 6GCL model, but not 
that of the isotropic six-state clock model. In the PD phase, large fluctuation appears, which is 
similar to the lOP on the 6GCL model. This large fluctuation made confusion in previous studies. 

The transition between the PD phase and the FR phase is of the first order. However, no 
anomaly has been found in experiments and in simulations. We have studied microscopic spin 
ordering in the chains along the c-axis to clarify the nature of the phase transition. That is, we 
study whether the change at the low temperature critical point is just a crossover from the PD 
type to the 2FR phase or a first order phase transition. We found a discontinuous change of the 
distribution of domain walls in the chains of the sublattices, and found the change is of the first 
order. However, the total number of domain wall changes very little at the point, which causes a 
very small singularity of the specific heat. 

From the viewpoint of the renormalization group,^^' the first order is also supported. Suppose 
the effective Hamiltonian of the present model is 



where $ is the complex field and $ is its conjecture. The renormalization group flow diagram of 
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this model is depicted in Fig.9. The 2FR phase and the PD phase correspond to Xq — oo and 
^ oo, respectively. Thus, the transition point between the PD and the 2FR phase is at Ae = 0. 
The nature of this transition at Ae = is similar to that of the XY model when the external 
field H change its sign. Below the transition temperature, the first order phase transition exists 
at ^ = 0. This first order transition of the XY model has a unique property that the correlation 
length diverges at ^ = 0. However, in the present STIAF model, due to the higher terms of 
^6n{n- > 2), it is expected that the transition becomes an ordinary first order transition where the 
correlation length is finite. Therefore, the transition between the PD phase and the 2FR phase is 
of the ordinary first order. 

We hope that the present study settles the long-studied problem of the low temperature tran- 
sition of the STAFI model. 
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Fig. 1. Schematic sublattice magnetization of (a) 2FR phase (b) 3FR phase (c) PD phase. 
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Fig. 2. Temperature dependence of the specific heat. 
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Fig. 3. The exchange interactions in the stacked triangular lattice. 
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Fig. 4. Order parameter space of the STAFI model. The solid line connecting a vertex and the center is one of the 
six low temperature states and dotted line between the neighboring solid lines is one of the six PD states. The 
three arrows denote the sublattice magnetizations: (mi, 7712,7713). The angles are assigned by the relation ?? 
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Fig. 5. Scahng plot of Mq around the low temperature phase transition pont {k^Tc/Jo ~ 0.188) for the various 
system size. 
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(a) (b) (c) 

Fig. 7. Clusters for the mean field approximations with mean fields. The gray circles represent the mean fields and 
white circles represent the spins of the cluster. 




Fig. 8. Temperature dependence of the sublattice magnetizations obtained by the mean field approximations with 
clusteres depicted in Fig. 7. 
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Fig. 9. Renormahzation flow diagram of the Ze symmetry model which project the t-Xe surface. NG and 3D-XY 
means Nambu-Goldstone fixed point and 3D-XY fixed point, respectively. 



